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Abstract: Z' gauge bosons arise in many particle physics models as mediators between 
the dark and visible sectors. We exploit dark matter complementarity and derive stringent 
and robust collider, direct and indirect constraints, as well as limits from the muon magnetic 
moment. We rule out almost the entire region of the parameter space that yields the right 
dark matter thermal relic abundance, using a generic parametrization of the Z'-fermion 
couplings normalized to the Standard Model Z-fermion couplings for dark matter masses 
in the 8 GeV-5 TeV range. We conclude that mediators lighter than 2.1 TeV are excluded 
regardless of the DM mass, and that depending on the Z'—fermion coupling strength much 
heavier masses are needed to reproduce the DM thermal relic abundance while avoiding 
existing limits. 
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1 Introduction 

One of the most exciting and tantalizing puzzles of modern cosmology and particle physics 
lies with the nature of dark matter (DM). DM comprises roughly 23% of the energy budget 
of the Universe and its existence has been so far inferred only from gravitational effects. 
The most compelling particle DM candidates are the so called Weakly Interacting Massive 
Particles (WIMPs), arising in a wide variety of well-motivated theories beyond the Standard 
Model (SM) such as supersymmetry and models with extra space dimensions. WIMPs can 
naturally account for the observed DM abundance in the framework of standard thermal 
freeze-out, and they are often within reach of current and future experiments (see Ref. [1-3] 
for extensive overviews on particle DM). 

Although the presence of DM has been ascertained only gravitationally, there are 
several observations which indicate that the direct detection of DM particles is around the 
corner. Direct detection experiments attempt to measure the recoil energy deposited by 
DM particles in underground experiments. The measured rate of scattering events, after 
the subtraction of background, singles out a preferred region in terms of the scattering 
cross section and DM mass (see Ref. [4] for a review of current status and prospects in this 
field). 

Indirect DM detection consists of the observation of DM annihilation products by 
satellite or ground-based telescopes. This is a promising and complementary avenue, since 
it extracts information concerning the DM distribution and annihilation cross section [5]. 
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Ref. [6] gives an account of the current status of indirect searches and of future observational 
avenues. 

Concerning collider searches, the typical signature of DM production is the presence of 
missing energy, since WIMPs simply escape the detector. Hence, collider searches for DM 
production typically comprise jets and missing energy and provide rather complementarity 
limits [7]. However, often times the most stringent limits do not come from missing en¬ 
ergy signatures, but instead from the production of SM particles with the mediator being 
produced in the s-channel. Looking for resonances in dileptons or dijets has proven to be 
a prominent way to constrain dark sectors. As we shall discuss further, those bounds are 
more stringent than the direct and indirect DM detection ones if the mediators interact 
with leptons. For recent leptophilic DM models see Refs. [8-10]. 

Several Dirac fermion DM models have been proposed, with DM particle masses in 
the GeV-TeV range, where the DM is coupled to the visible sector via Z' gauge bosons, 
arising from a U{l)x gauge sector. Such gauge bosons have been extensively studied in 
many different contexts [7, 8, 11], but our focus here is on the dark matter phenomenology. 

Instead of choosing a particular model, we perform here a comprehensive, model- 
independent analysis. Aside from an overall rescaling, we assume that the Z' has similar 
couplings to SM Z, an assumption studied in the previous literature [7, 11, 12]. By studying 
regimes in which the Z'-SM fermion couplings are suppressed by an overall factor, we 
implicitly encompass several SM extensions such as sequential Z' models [13], 3-3-1 models 
[14, 15], E-6 models [16], Left-Right [17], etc. We present results for a broad range of 
particle DM masses, specifically: 8 GeV, 15 GeV, 50 GeV, 500 GeV, 1 TeV, and 5 TeV, 
with Z' — X — X couplings ranging from 10“^ — 10. 

Due to the presence of a spin 1 mediator that couples to the muon, important con¬ 
straints potentially arise from contributions to the muon anomalous magnetic moment. We 
use here a new public code for computing the muon magnetic moment to derive limits on 
the mediator mass described in Ref. [18]. Since only the muon’s interactions with the Z' 
enter into this calculation, the muon magnetic moment offers quite interesting bounds in 
the Z'-DM suppressed couplings regime where direct and indirect detection experiments 
lose sensitivity. 

As mentioned, extensive literature exists on the subject topic we focus on here. Our 
work is similar in vein to that presented in Ref. [7-10, 14, 19-24]. Ref. [7] focused on the 
leptophobic regime, whereas Refs. [8-10] discussed the leptophilic scenario. In Refs. [14, 
19, 20] the DM phenomenology was studied in a particular model, and Ref. [21] focused 
on resonances scenarios. Ref. [22] performed a similar analysis to ours with an older data, 
without computing indirect DM detection and g-2 bounds, and assumed Z' — up quarks 
and Z' — down quarks to be the same, whereas in Ref. [23] the kinetic mixing window has 
been explored. In Ref. [24] the authors discussed the very light mediator regime only. 

The present study, however, differs from the mentioned previous studies in several key 
ways: 


• We use a generic parametrization to describe the Dirac fermion DM - SM fermion 
interactions in the context of the Z' portal. 
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• We derive up to date spin-independent and spin-dependent limits using LUX and 
XENONIOO data. 

• We derive indirect detection limits, including the important annihilation channel to 
on-shell Z's when mz’ < 

• We carry out a detailed and comprehensive collider analysis, using up to date dilepton 
limits. 

• We obtain limits from the muon magnetic moment on the Z' mass. 

We begin our study with an outline of the Z' DM portal we focus on. 

2 Z' Portal Dark Matter Model 

The existence of new gauge symmetries is among the best motivated extensions to the SM. 
In particular, in the low energy limit of many new-physics scenarios, an additional heavy 
electrically neutral gauge boson, Z', generically possesses sizable couplings to SM quarks 
and leptons [25]. DM that interacts with the SM purely through electroweak interactions 
is strongly constrained [26], and therefore a dark sector that communicates with the SM 
via a new Z' remains a very natural scenario. In this paper, we study a generic Z' model, 
and parametrize the Z'-fermion couplings so that the Z' couples similarly to the SM Z 
boson (see e.g. Ref. [11]). This regime is also known as a sequential Z'. In this case, the 
crucial differences between the Z' and the Z are their masses and the Z' couplings to new 
particle species (DM for example). Additionally, we add a universal scale factor “a” in such 
a way that when a = 1, then the Z' couples to the SM fermions equivalently to the SM Z, 
whereas when a differs from 1 we are accounting for suppressed couplings. In addition to 
the unsuppressed case (a = 1), we analyze the case where a = 0.5. In other words, the Z' 
couplings to the SM are universally scaled down (relative to the purely sequential case) by 
a factor of 50%. 

There are many types of DM that can interact with Z' bosons. However, to simplify the 
discussion we restrict ourselves to DM composed of a single Dirac fermion, y. Obviously, 
the introduction of a Z' boson requires, in principle, the inclusion of additional fermions 
to cancel the triangle anomalies. Since we focus on the collider and DM phenomenology of 
the Z' portal we will remain agnostic about the specific scenario for anomaly cancellation. 
For a review of this topic concerning Z' models see Ref. [27]. Then, a model-independent 
approach involves parametrizing the simple tree-level Z' interactions as 

CD Z'^ [X7^ (s-^v + X + ax (^/v + Qfal^) /] , (2.1) 

/eSM 

where the sum is over all quarks and leptons of the SM (including neutrinos). We note that 
for Majorana DM, the most significant change would be when considering constraints from 
direct detection since then = 0. For Dirac DM, these vector interactions are strongly 
constrained by current measurements of spin-independent scattering with nuclei [28] , as will 
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Figure 1. a,b) annihilation channels; c) scattering diagram; d) muon magnetic moment diagram. 


be discussed in Sec. 3. To further simplify the parameter space, we will assume x’s vector 
and axial interactions with the Z' are of equal magnitude and in particular that = 9xa- 
For this reason, we will often refer to the Z'—DM coupling simply as = fl'xa- 

Choosing instead g-^^ = —gxa would not significantly change our results. Moreover, if one 
departs from g^.^ ~ g^a couplings assumption, in particular in the scenario which vector 
couplings are suppressed along with direct detection limits, the collider bounds we derive 
further will remain at the TeV scale, keeping our overall conclusions unchanged. 

In the case of a purely sequential Z' (a = 1), the couplings of Z' to the SM fermions 
are exactly the same as the SM Z. Even though this type of Z' is only expected from gauge 
theories where the Z' and the SM Z have different couplings to new exotic fermions, it is 
a useful reference model and is analogous to new states of the SM Z in theories of extra 
dimensions at the weak scale. However, we will additionally consider the parameter space 
where the couplings gfy, gja are all universally suppressed by some constant factor “a” 
introduced in Eq. 2.1 [12]. More precisely, for the Z' couplings to SM quarks and leptons, 
we set 

— g 

gua = -j- (tan + cot ffiu) 

g 

gda = ^ (tan 01 ^ + cot6»^) 

g 

9ia = ^ {tan6w + cot 9ui) 

— g 

gua — ^ (tan Oyj -j- cot , (2.2) 

where corresponds to up-type quarks, down-type quarks, charged leptons, and 

neutrinos, respectively. Furthermore, e = gy,sin6y, is the electromagnetic coupling, 0^ is 
the Weinberg angle, and a = 1 corresponds to an exactly sequential Zb Our results will 
focus on the values a = 1 and 0.5. 

We emphasize that our model independent approach maps onto several models dark 
matter models such as, but limited to, [13-17] 

3 Direct Detection 

Given the simplihed model of Eq. (2.1), y may scatter off of quarks inside nuclei through 
the t-channel exchange of a Zb Eurthermore, if the Z'-y interactions are sufficiently large, 


—e /5 

5«v = ^ ( 2 tan 9yj - cot 6'^, 

e A « 

9dv = ^ 1 2 tan 9y, - cot 6 

g 

9iv = I (3 tan 6*^ - cot 0^) 

g 

guv = ^ (tan6»^„ -7 cot 0^) 
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then these scattering rates are subject to current bounds from direct detection experiments 
like LUX [28] and XENONIOO [29]. In this section, we will give the approximate form for 
these rates. In the limit of small momentum transfer, Z' exchange can generically lead to 
spin-independent (SI) scattering with nucleons n. 


a 


SI (per nucleon) 


vr 


■2^/prot + (^ ~ Z)fneut 


A 


/prot — 2 {“^duv + 9dv) 


nir 


/nent = ^ (5 uv “t“ ‘^Qdv 


rrir 


(3.1) 


and spin-dependent (SD) scattering, 


^SD (per neutron) 


3n /^^neut 




TT 


m 


z' 


(Ar* + 


1 2 


(3.2) 


where is the WIMP-nucleon reduced mass, Z, A are the atomic number, atomic mass 
of the target nucleus, respectively, and are the quark spin fractions of the neutron. 

We will take these to be = —0.42, A))®"* = 0.85, A“"* = —0.08 [30] and require 

the spin-independent and spin-dependent rates to be below current published limits from 
LUX [28] and XENONIOO [29], respectively. Note that for spin-dependent rates, we con¬ 
sider scattering with neutrons since these limits are usually stronger. 


4 Indirect Detection 


If X is thermally populated through the interactions present in Eq. (2.1), it is possible 
that its annihilation rate into SM fermions is still large today. In fact, annihilations into 
final state quarks and leptons can lead to a diffuse emission of high energy gamma-rays 
through the processes of neutral pion production and final state radiation, respectively [3]. 
Indirect detection measurements, such as those performed by the Eermi satellite, seek to 
observe these emissions in the galactic center (GC) of the Milky Way. Of course, certain 
astrophysical assumptions must be made (e.g. the nature of the DM density profile), 
but even conservative ones can lead to significant constraining power in the annihilation 
rate today, {av) [31]. In this section, we briefly describe the analytic forms for these 
annihilations and the method by which we use the results of Fermi observations to constrain 
the parameter space of our Z' models. 

From the Lagrangian of Eq. (2.1), the non-relativistic form for the annihilation cross 
section into a pair of SM fermions, /, through an s-channel Z' (assuming that ^ i^f) 
is 


o-u {xx ff) 


a‘^nc\ll- 
27rm^, (4m^ — 

I 2 2 4 /o 2 I 2^ 

+ 9xv9fv'^z' -|- rrij-) 


9fa 


‘^9lym%, {m\ - rrij) + 9ja'fnj (4m: 


ml,f 

(4.1) 
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where v is the relative velocity of the annihilating DM pair and ric is the number of colors 
of the final state SM fermion. When near resonance, the form of Eq. (4.1) must be properly 
modified by including the Z' width, E^/, which from Eq. (2.1) has the form 


Tz' = ^ [mz' - ^mf ) 
/eSM 


+ 6 {mz' — 2m^) 


mz' 

127r 



1 - 


4m^ 


nir 


(4.2) 


where 9 is the unit step function. In Eigs. 4-9 we present a light blue excluded region that 
represents the violation of the perturbative limit Vz' jvaz' <0.5. Obviously, in this region 
the narrow width approximation fails [32], and our limits are not applicable. We show such 
regions for completeness though. Experimental searches for spin 1 bosons at the LHC use 
as benchmarks the Sequential Z' Model and other similar models, which predict a narrow 
width for the Z' boson. In other words, those searches assume that a narrow resonance can 
be described by a Breit-Wigner line-shape when studying the distribution of the invariant 
mass of the decay products, and this assumption definitely breaks down along with the 
collider bounds for E^z/M^/ >0.5 [32]. Since the width of the Z' will largely be controlled 
by decays into DM pairs, large couplings can still be present if Mz' < 2m^. Similarly, the 
relic abundance calculation also requires a narrow width approximation. Hence, hereafter 
we will show the excluded large width region in blue in Eigs. 4-9. 

The variety of channels increases further if annihilation directly into pairs of on-shell 
Z' bosons becomes kinematically accessible. In this case, the non-relativistic form for the 
annihilation cross section (not including the decay of the Z'') is 


av [xx ^ Z'Z') 




1 - 




mf 


3/2 


1 - 


mr 


-2 


2m |) 






+ (ffxv + 9 


XO- 


^9-xy9-)^a 


) 


(4.3) 


Once produced in this manner, the Z' bosons will decay to pairs of SM fermions. Eq. (4.3) 
must then be modified by including the appropriate factors of branching fractions for 
Z' —^ //. Therefore, when > mz', these 2-stage, or cascade, annihilations can lead 
to gamma-rays. Compared to the direct annihilation case, the Lorentz boost between the 
Z' and SM fermion rest frames generally leads to a widening in the energy range of the 
gamma-ray spectra for cascade annihilations. 

Eor direct annihilations, as in Eq. (4.1), we directly utilize PPPC4DMID to calculate 
the spectrum of gamma-rays [33]. However, when considering the cascade annihilations of 
Eq. (4.3), the direct spectra obtained by PPPC4DMID must be properly integrated over a 
finite energy range, corresponding to the range of boosts that the final state fermions are 
kinematically allowed to obtain [34]. 

Given certain astrophysical assumptions, indirect detection measurements can only 
constrain the quantity •: where {av)j and N-^j are the annihilation rate 

and number of photons emitted (with energy in some specified range) per annihilation 
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for some final state /, and the factor of | takes into account that Dirac DM is not self¬ 
conjugate. In particular, we use the constraints obtained in Table I of [31], which gives 
95% CL upper bounds on this quantity for different assumptions of the DM profile and 
gamma-ray energy bin.^ In comparing to our theoretical prediction for ^ Ylf /-^ 7 ,/) 
we utilize bounds that are derived using a Navarro-Frenk-White (NFW) profile, and in 
each scenario that we will consider, we use the gamma-ray energy bin that gives the most 
stringent upper limit in the parameter space of interest. 

In the two previous sections, the direct and indirect detection observables were de¬ 
scribed. We now briefly review and compute the correction to the muon magnetic moment 
arising from our model. 


5 The Muon Anomalous Magnetic Moment 


The muon magnetic moment is one of the most well measured observables in particle 
physics. There is a long standing discrepancy between theory and experiment of about 3.6a. 
The large uncertainty surrounding the theoretical calculation blurs the excess. The current 
deviation is Aa^ = 295T81 x 10“^^. Out of this ±81 x 10“^^ error, ±51 x 10“^^ is theoretical, 
±39 X 10“^^ rising from the lowest-order hadronic contribution and ±26 x 10“^^ from the 
hadronic light-by-light correction. Since the Z' interacts with the muon, it also gives rise 
to corrections to the muon magnetic moment through Fig. 1(d), with a contribution found 
to be [18, 35], 


A alvPAx) + g^aPaix) 

Jq "^^(1 - a:)(l - X^x) ± A2x’ 


(5.1) 


where A = m^/Mz' and 


Pv{x) = 2x^(1 — x) 

Paix) = 2x(l — x) • (x — 4) — 4A^ • x^. (5.2) 


In the limit Mz' 3> this integral can be simplified to, 


Aa^(Z') 


ml 

47r2M|, 


^ 2 
^9fiv 



(5.3) 


which agrees with Ref. [35]. Notice that, depending on the relative values of the vector 
and axial couplings, the contribution can be either positive or negative. In this work, 
the contribution is always negative, since the Z' has similar vector and axial couplings to 
the SM Z. Thus, we can place the Icr bounds based on the aforementioned error bar. In 
Figs. 3-8, when we refer to the g-2 bound, we will be forcing the Z' contribution to be 
within the the la error bar. 

We now will turn our attention to the collider phenomenology, focusing mostly on the 
dilepton searches. 

^In comparing Table I to Fig. 10 of [31], we find a discrepancy, which results in the need to rescale the 
values in Table I by a value of approximately 2.8 (for an NFW profile). This has been confirmed with the 
authors of [31]. 
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6 Collider Bounds 


6.1 Bounds from the search for new resonances in dilepton events 

Several extensions of the SM predict the existence of new gauge bosons as a heavy Z', 
which can be detected by looking for a resonance in dijet or dilepton invariant masses. The 
most stringent constraints to date on heavy Z' bosons decaying to charged leptons come 
from the 8 TeV LHC data [36]. 

We find that, in fact, constraints from dijet and monojet searches for a Z' are compet¬ 
itive only in the leptophobic scenario studied in [12] confirming those results. As long as 
the couplings of the Z' to charged leptons are non-negligible, the process pp —)• Z' —>■ 

£ = e, //, becomes more relevant compared to dijets and mono-A processes. 

Concerning the comparison between dijet and dilepton processes, as the couplings 
to leptons and quarks are of electroweak order, the production rates are not drastically 
different, with quarks being favored due they color multiplicity. However, the backgrounds 
to dijets involve QCD pair production and are much larger than SM Drell-Yan background 
in the case of leptons. Hence, monojet searches can be competitive only for large Z' 
couplings to neutrinos or to DM. However, a scenario with large couplings to neutrinos 
and small couplings to charged leptons is not as likely. Furthermore, a very large coupling 
to DM is not favored by direct detection data. The effect due to the mass and couplings 
of the DM particle is indirect concerning the way a dilepton search might constrain a Z' 
DM model, but not to the way a specific model might fix the couplings between Z' and 
the SM fermions. In the first row of Fig. 2 we show the branching ratio of Z' into jets, 
charged leptons and an invisible mode comprising dark matter decays and neutrino decays. 
The branching ratios are constant as a function of the Z' mass except near the threshold 
to the dark matter channel. Note that the invisible mode is always sizable because of the 
neutrinos. 

The situation is similar as far as the variation of the dark matter mass is concerned, 
as shown in the second row of Fig. 2. The branching ratios of all modes are constant away 
from the threshold region and the relative size of the branching ratios are controlled by 
the couplings a and g-^. To understand how the couplings affect the branching fractions, 
we display in the third and the fourth rows the variations against and a, respectively. 
We kept the Z' and the dark matter masses fixed at a representative point away from any 
threshold region: = 2 TeV and = 130 GeV. It is clear from these plots that the 

fraction of decays into dark matter and SM fermions are determined by the size of the 
couplings we just mentioned. 

To summarize, for a hxed Z' mass, the branching ratios into charged leptons BR{Z' —)• 
and into jets BR{Z' —?■ jj) decrease as soon as the DM channel becomes available. 
Also, decreasing the coupling g^. between the DM and Z' increases the branching faction 
into leptons and quarks, the same effect occur when we increase the a coupling. However, 
the branching ratio to charged leptons are much less sensitive to variations of the models 
parameters compared to the jets and dark matter channels, varying slightly around the 
10% level as can be seen in all panels of Fig. 2. 
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Figure 2. Branching ratios of Z' into jets (dotted lines), leptons (solid lines) and an invisible 
(neutrinos plus DM) mode (dashed lines) as a function of the new gauge boson mass in the first 
row, the DM mass in the second row, the Z'-DM coupling in the third row, and the scaling 
factor a of the coupling between Z' and the SM fermions in the fourth row. 


The branching ratios to jets are considerably smaller than those found in the leptopho- 
bic scenario of [12] in the entire Z' mass range for all dark matter masses. The branching 
ratio to an invisible final state is enhanced compared to a leptophobic scenario due to the 
appearance of a neutrino channel. As we see in Fig. 2, the branching ratio to DM plus 
neutrinos dominate at some points of the parameter space when ;§> a. This is almost 
twice the typical branching ratio to an invisible channel we might expect in the leptophobic 
case when ~ 1 [12]. Yet, the impact on the monojet channel is mild. The very hard cut 
on missing energy imposed by the LHC studies [37] selects mainly events with DM instead 
of neutrinos. 

In summary, allowing sizable Z' couplings to charged leptons softens the constraints 
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found in the leptophobic case studied in [12]. In that work, for g^<l,Z' masses up to ~ 2 
TeV for light DM could be discarded based on dijet and monojet searches at the Tevatron 
and LHC, but those limits get much less stringent as decreases. As we are going to 
show, the limits from searches for a dilepton resonance do not soften in the small Z'-DM 
coupling regime. As we now show, dilepton searches are more efficient to exclude regions 
of the parameter space of dark Z' models which present sizable lepton couplings. 

Dileptons at the 8 TeV LHC - The ATLAS search 


To evaluate the impact of the 8 TeV LHC search on dilepton resonances after 20.3 fb~^ 
of integrated luminosity for the electron sample and 20.5 fb~^ for the muon sample [36], 
we simulate the process 

pp^ Z' , £ = e,g (6.1) 

plus up to two extra jets using MadGraphS [38]-FeynRules [39], clustering and hadronizing 
jets with Pythia [40], and simulating detector effects with DelphesS [41]. Soft and collinear 
jets from QCD radiation generated by Pythia are consistently merged with the hard radi¬ 
ation calculated from matrix elements in MLM scheme [42] at appropriate matching scales. 
We adopt the CTEQ6L parton distribution functions computed at fip = PR = Mz'- 

The signal events were then selected with the same criteria adopted in [36] shown 
below 


Priei) > 40 GeV , pt{& 2 ) > 30 GeV , [rye] < 2.47 (6-2) 

Pt{pi) > 25 GeV , pt{P 2 ) > 25 GeV , [ry^] < 2.47 (6.3) 

128 GeV < Mu < 4000 GeV (6.4) 

Where is the hardest (second hardest) lepton in the event, and Mu the invariant 

mass of the lepton pair. The signal acceptance times efficiency found in our simulations 
are similar to those presented in [36]. 

All the backgrounds simulations to pp ^ Z' ^ l'^ were taken from [36]. To place 

limits on a Z' model we calculated a statistic at the 95% confidence level based on Mu 
measured in [36] in 6 invariant mass bins: 110 — 200 GeV, 200 — 400 GeV, 400 — 800 GeV, 
800 — 1200 GeV, 1200 — 3000 GeV, and 3000 — 4500 GeV. We show in Fig. 3 the number of 
events for the assumed luminosity in each invariant mass bin for the total background 
and signal for various values of and the production cross section as a function of M^. 

The limits obtained from the dilepton invariant mass depend weakly on the DM mass 
and softens as g^^ approaches 1. This can be easily understood reminding that the branching 
ratio to charged leptons are only moderately dependent upon the relevant parameters M^, 
g-^ and a as can be seen in Fig. 2, except in the strong g^ regime. Overall, the limits 
obtained on the Z' mass are close to those found in the experimental study of the ATLAS 
collaboration [36]. 

We do not take systematic uncertainties into account in the computation of the limits 
on the Z' mass. For this reason, the colliders bounds are not likely to be conservative. 
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Figure 3. The dilepton invariant mass distribution for the total background (dashed line) in six 
bins taken from [36] and the signal for four different Z' masses: 500 GeV, 1, 2 and 3 TeV. In the 
inset plot we show the total cross section in pb as function of Mz'- The integrated luminosity 
assumed in this plot is the same of the experimental study [36]. 

However, the effects from mis-modelling the shape of the signal and background distribu¬ 
tions are not expected to be large due the coarse binning used in the experimental study, 
that is it, it is not necessary a large number of free parameters to fit the distributions. 
For example, in [43] a detailed analysis involving a bin-by-bin systematic uncertainty is 
shown to have only a moderate effect in the collider bounds even for a finely binning Mu 
distribution. 

7 Dark Matter Complementarity 

We have discussed a collection of DM constraints stemming from direct, indirect and col¬ 
lider searches and the muon magnetic moment that are relevant for the Z' portal scenario 
under consideration. In this section, we combine these limits and outline the viable ver¬ 
sus excluded region of the parameter space. Also, we explicitly show how critical DM 
complementarity is. We focus on 8 GeV, 15 GeV, 50 GeV, 500 GeV, 1 TeV, and 5 TeV 
WIMP masses, motivated by the recent tentative DM signals and the available LHC data 
[44]. The direct detection constraints are derived using LUX 2013 and XENONIOO 2012 
results, whereas the indirect detection limits are obtained using the Fermi Galactic Center 
(Fermi GC) limits from [31]. Moreover, the collider limits from dilepton resonance searches 
are derived using 8 TeV ATLAS data as discussed above. Additionally, we included the 
perturbative limit on the Z' width, which removes the region with Vz' jMz' > 0.5 [32] as 
discussed in the previous section. Lastly, the muon magnetic moment limits are obtained 
using the public code [18], where we have simply adapted it to our model. We emphasize 
that the collider, direct and indirect detection bounds are placed with 95% C.L, whereas 
the muon magnetic moment is a \a limit. 

In Figs. 4-9 we exhibit the parameter space ruled out by those combined limits, where 
the importance of DM complementarity is made clear for DM couplings ( 5 ^) larger than 
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Figure 4. Blue horizontal line is LHC exclusion. Everything below the curve is excluded. Gray 
horizontal line is the Icr bound from the g-2. In red, green and pink we show the LUX spin- 
independent, Fermi Galactic Center, and XENON spin-dependent limits respectively. The black 
curve yields the right abundance. The light blue region region corresponds to where the perturbative 
limit on the Z' width is violated. Left: 8 GeV WIMP with Z' = Z■,Right: 8 GeV WIMP with a 
Z'-fermions couplings 50% suppressed compared with the SM Z. The somewhat wavy aspect of the 
limit for large couplings occurs due to the coarse binning of the dileptons invariant mass distribution 
delivered by the experimental study. 



Figure 5. We have used the same color pattern throughout. Left: 15 GeV WIMP with Z' = Z\ 
Right: 15 GeV WIMP with a Z'-fermions couplings 50% suppressed compared with the SM Z. 


unity. 

In the left panel of Fig. 4, we can clearly see even without the inclusion of dilepton 
bounds that 8 GeV WIMPs are completely excluded by the combined direct detection and 
the muon magnetic moment limits. Since the Z' is a Z-like gauge boson, the correction to g- 
2 is negative and therefore a limit of M'^ > 135 GeV is found using Eq. (5.3). As predicted, 
the muon magnetic moment places complementary limits in the regime of suppressed Z' — 
X ~ X couplings, where direct and indirect detection ones become loose. In addition, 
dilepton data implies > 2.5 TeV for g^. < 0.2. For couplings larger than unit, the 
collider bounds are strongly weaken. 

As gy- increases, the branching ratio to leptons drops and the collider bounds from 
dileptons are evaded for all Z' masses. In the region of g^^ where those bounds weakens 
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Figure 6. Left: 50 GeV WIMP with Z' = Z;Right: 50 GeV WIMP with a Z'-fermions couplings 
50% suppressed compared with the SM Z. 




Figure 7. Left: 500 GeV WIMP with Z' = Z;Right: 500 GeV WIMP with a Z'-fermions couplings 
50% suppressed compared with the SM Z. 



Figure 8. Left: 1 TeV WIMP with Z' = Z]Right: 1 TeV WIMP with a Z'-fermions couplings 50% 
suppressed compared with the SM Z. 
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Figure 9. Left: 5 TeV WIMP with Z' = Z;Right: 5 TeV WIMP with a Z'-fermions couplings 50% 
suppressed compared with the SM Z. 


the shape of the curves seen in Fig.4, are somewhat wavy due the coarse binning of the 
dileptons invariant mass distribution delivered by the experimental study. When the Z' 
mass suddenly populates another bin there’s a fast weakening of the bounds causing a 
tooth-like curve. This feature is repeated for all DM masses. 

In the right panel, we have universally suppressed the Z'-fermion couplings by 50%, i.e 
we set a = 0.5 in Eq. (2.1). This suppression does not help much, since the entire region of 
the parameter space that sets the right abundance is ruled out, but the complementarity 
between DM and collider searches starts to become relevant for large couplings. In this 
case g-2 demands > 67 GeV, and dileptons data requires >2.1 TeV. In this regime 
is clear the importance of exploiting complementary searches. 

In Fig. 5, it is already quite noticeable the complementarity among direct, indirect, 
collider searches, and g-2 mostly for WIMP masses of 15 GeV as we are now away from 
the energy threshold in direct detection experiments. Again, the dileptons constraints are 
the most stringent ones, ruling out Z' masses below 2.5 TeV (2.1 TeV) in the left (right) 
panel for couplings smaller than 0.1. For larger couplings the direct detection limits are 
the most stringent ones. The bounds from the muon magnetic moment will follow the one 
in Fig. 4 since it is independent of the DM mass. 

Similar conclusions are found for = 50 and 500 GeV in Figs. 6-7. In Fig. 8, for 
rriy- = 1 TeV, we start seeing non-collider bounds becoming more competitive for large 
DM couplings and some region of the parameter space that reproduces the correct relic 
abundance surviving all constraints. Interestingly, the combination of g-2, direct, indirect 
and collider data excludes < 10“^ in the left panel and < 10“^ in the right panel 
and favors only heavy mediators. Lastly in Fig. 9 for = 5 TeV, we see that this feature 
from Fig. 8 continues and now a much larger region obeys all limits. Albeit, interestingly 
the perturbative limit on the Z' width removes the region g^ ~ 10 and Mz' ~ 10 TeV. 

In summary, when one properly takes into account dilepton, g-2, direct and indirect 
detection data, a Z' lighter than 2.1 TeV is ruled out, independent of the DM mass. One 
therefore needs vector mediators heavier than 2.1 TeV to accommodate a Dirac fermion 
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DM particle. We emphasize that depending on the Z'-fermion coupling strength masses 
much heavier than 2.1 TeV might be required, regime As we increase the DM mass, the 
indirect detection limits become more relevant . In Figs.8-9 the indirect detection limits are 
some times stronger than the direct detection ones. Moreover, the dilepton limits which are 
usually the most relevant ones, become secondary as larger Z' — DM couplings are used, 
showing a strong degree of complementarity between those searches. We emphasize that we 
have used a Z'-fermion parametrization which is similar to the SM Z one. Since we studied 
different DM masses,Z'-fermion couplings and a large range for the Z'-dark matter coupling 
magnitude we have covered several U{l)x extensions of the SM, in a model independent 
fashion. It would be possible to alleviate those constraints by advocating the presence of 
a leptophobic Z' as in [12], or a Majorana DM fermion such as in Ref. [45], or possibly a 
pure axial Z'-fermion interactions such as in Ref. [46]. 

If one departs from the ~ g^a couplings assumption, specially in the regime that 
vector couplings are dwindled, the direct detection limits are automatically suppressed, 
whereas the collider bounds will remain in the TeV scale, keeping our overall conclusions 
unchanged. We point out that the DM phenomenology in this suppressed scenario has 
already been studied in Ref. [11]. 

We point out that if one had used a different parametrization scheme for the Z'- 
fermion interactions mild changes are expected and the general statement that Z' portal 
only allows heavy mediators is still valid. Obviously, a key assumption made throughout 
this work pertains to the strength of the Z'-fermion coupling. Since we have normalized 
our results in terms of the SM Z coupling strength, one could evade those limits advocating 
much more suppressed couplings. On the other hand, one could also use larger Z'-fermion 
couplings as long as the Z' is sufficiently heavy in order to avoid the discussed limits. 

8 Conclusions 

Motivated by potentially exciting direct and indirect detection signals, in this work we ex¬ 
ploited DM complementarity in the context of the Z' portal using a generic parametrization 
of the Z'-fermion couplings with Dirac DM. We performed a detailed analysis of collider, 
direct and indirect DM detection data as well as the muon magnetic moment to outline 
the viable vs. excluded region of the parameter space. 

A high degree of complementarity is observed at several different levels: The muon 
magnetic moment provides complementary limits in the regime of suppressed x — x — Z' 
couplings, since then only the Z' — g coupling strength is relevant. Indirect detection 
limits, on the other hand, trace back the region of the parameter space that sets the 
right abundance, regardless of the coupling strength, and rules out WIMP masses below 
15 GeV (except for some viable parameter space very close to resonance). Since we are 
using a generic parametrization of the Z'-fermion interactions, both spin-independent and 
spin-dependent scattering exists in this model, but we have shown in Figs. 4-9 that the 
LUX spin-independent limits overwhelms the spin-dependent ones due to the well known 

enhancement from coherent nucleon scattering. Furthermore, we have shown that 
for Z' — DM couplings smaller than unit, the dilepton bounds are the most stringent 
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constraints derived in this work, whereas for larger couplings the indirect and most often 
the direct detection ones are the most restrictive, emphasizing the importance of DM 
complementarity. They almost exclude the entire region of the parameter space of the 
model that sets the right abundance allowing only heavy mediators attached to heavy DM. 
We emphasize that the light DM regime is ruled not matter how heavy the Z' mass and 
that our limits would be obviously stronger if larger Z'-fermion couplings were used. 

In conclusion, after varying the DM mass from 10 GeV up to 5 TeV, and after consider¬ 
ing suppressed couplings, we have demonstrated that Z' mediators lighter than ~ 2.1 TeV 
are excluded regardless of the DM particle, and, depending on the Z' — fermion coupling 
strength, much heavier masses are needed to reproduce the DM thermal relic abundance 
while avoiding existing limits. 
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